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^ ; Abstract 
> ■ 

[-^ ■ We have calculated the radiative corrections to the Dalitz plot of baryon semileptonic decays with angular 

T-H ! 

Q^ ■ correlation between polarized emitted baryons and charged leptons. This work covers both charged and 

^. 

r I , neutral decaying baryons and is restricted to the so-called three-body region of the Dalitz plot. Also it 

o 

QQ I is specialized at the center-of-mass frame of the emitted baryon. We have considered terms up to order 



{a/iT){q/Mi)^, where q is the momentum transfer and Mi is the mass of the decaying baryon, and neglected 



^ , terms of order {a/'ii){q/MiY for n > 1. The expressions displayed are ready to obtain numerical results, 

H ■ 

Cu \ suitable for model-independent experimental analyses. 

PACS numbers; 14.20.Lq, 13.30.Ce, 13.40.Ks 



I. INTRODUCTION 

Currently, experiments on spin 1/2-baryon semileptonic decays (BSD), A — > Biui, where the 
polarization §2 of the emitted baryon B is observed are underway yj]. The analysis of these 
experiments requires the inclusion of radiative corrections (RC) to the Dalitz plot when §2 is 
nonzero. Our previous work ||2|] does not cover this case. It is the purpose of this paper to produce 
such RC. 

There are several requirements that must be met. In order to keep experimental analyses model 
independent it is necessary that RC are model independent themselves. There are many possible 
charge assignments to A and B and RC should be calculated so as to cover all the expected 
assignments. The charged lepton £ should be allowed to be an e^, /x^, and even r^ as the case 
may be. Since RC depend on the form factors present in the uncorrected decay amplitude, it 
is also necessary that they be cast into a form that can produce numerical results which are not 
compromised by fixing the form factors at prescribed values. 

The model independence of RC is achieved by following the generalization for hyperons [|30 of 
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the treatment of virtual RC in neutron beta decay [140 and of the application of the Low theorem 
151] to the bremsstrahlung RC developed in Ref. [16[]. There are six different charge assignments 
predicted by the light and heavy quark content of A and B. To cover all these cases it is necessary 
to know only the RC to the neutral decaying baryon (NDB) A^ -^ B^i'ui and to the charged 
decaying baryon (CDB) A^ -^ B^i^Ui cases. The other possibilities are obtained using the RC 
of the latter two [170. The cases £ = e^ , /x^, r"*^ are included in the RC by keeping the mass m 
of i uncompromised all along the calculation. In order to produce numerical values of RC that 
are practical to use in the Monte Carlo simulation of an experimental analysis and that are not 
committed to fixed values of the form factors of the weak vertex, one can numerically calculate 
the RC to the coefficients of the quadratic products of form factors that appear in the theoretical 
differential decay rate of the decay being measured. 

Since current experiments are medium-statistics (of the order of thousands of events) experi- 
ments and in order to keep the effort of calculating RC within convenient bounds, we shall consider 
contributions of order («/7r)(g/Mi)", with n = only and neglect orders with n = I and higher. 
Here q is the four-momentum transfer and Mi is the mass of A. Also we shall exhibit our results in 
a form where the integration over the real photon variables are ready to be performed numerically, 
except for the finite terms that accompany the infrared divergence of the bremsstrahlung RC which 



will be given analytically. The virtual RC will be given fully analytically. Our final result will be 
specialized to the center-of-mass frame of the emitted baryon B. 

In Sec. II we introduce our notation and conventions and discuss in detail the boundaries of the 
Dalitz plot in the center-of-mass frame of B. We shall specialize our calculation to the three-body 
region of this plot. Section III is devoted to the model-independent calculation of virtual RC. We 
will see that they can be put formally in the same form of our previous work, although now they 
will be functions of the energies E of i and Ei of A in the center-of-mass frame of B. The rather 
long expressions containing the form factors that appear in these corrections are exhibited in full 
in Appendix A. The bremsstrahlung RC are obtained in Sec. IV also in a model-independent form. 
However, the detailed discussion of its infrared divergence and the finite terms that accompany it 
is presented in Appendix B. In Sec. V we collect our results in a final form and we discuss their 
numerical use. We will cover the NDB and the CDB cases but we will exhibit only the calculation 
of the CDB case and limit ourselves to present the final results for the NDB case. Section VI is 
devoted to a brief discussion of our results. 

II. DALITZ PLOT IN THE CENTER-OF-MASS FRAME OF THE EMITTED BARYON 

For definiteness, let us consider the BSD 

A- ^B° + e~ + ue. (1) 

The four-momenta and masses of the A^, B^, i^, and Vi, will be denoted by pi = (_E'i,pi), 
P2 = (-E'2,P2), I = {E,\), and p,y = (E°,p,,), and by Mi, M2, m, and m^, respectively. The 
reference system we shall use is the center-of-mass frame of B. Accordingly, E2 = M2 and 
P2 = 0. It must be kept in mind that all other variables are referenced to this frame now. There 
should not arise any confusion with our previous work. A vanishing neutrino mass mj, will be 
assumed. Additionally, the direction of a vector p will be denoted by a unit vector p and whenever 
the expressions involved are not manifestly covariant, quantities like pi, I, or p^, will also denote 
the magnitudes of the corresponding three-momenta, unless stated otherwise. 

The uncorrected transition amplitude Mq for process ([U) is given by the product of the matrix 
elements of the baryonic and leptonic currents, namely, 

C 

Mo = -^[uBiP2)Wf,{pi,p2)uAipi)][ue{l)Of,v„{p^)], (2) 



where ua, ub, ue, and v,^ are the Dirac spinors of the corresponding particles and W^ is the weak 
interaction vertex given by 



W,{puP2) = /i(g^)7. + /2(g^)a,.^ + /3(g^)^ 



+ 



Ml ^^^^ 'Ml 



9i{(f)i^i + g2{q^)(yi,uT7- + 9z{q 



75- (3) 



Here O^ = 7^(1 + 75), q = pi — p2 is the four-momentum transfer, and fi{q'^) and gi{q'^) are the 
conventional weak vector and axial- vector form factors, respectively, which are assumed to be real 
in this work. In Eq. ^ we have omitted the Cabibbo-Kobayashi-Maskawa factors. They should 
be inserted once decay ^ is particularized. 

To cover the observation of the polarization of 5, its spinor is modified through the replacement 

ub{P2) -^ T.{s2)ub{p2), (4) 

where S(s2), the spin projection operator, is given by 

V/ N I-75/2 .^. 

S(s2) = , (5) 

and the polarization four- vector S2 satisfies the relations S2 ■ S2 = —1 and S2 ■ P2 = 0. In the 
center-of-mass frame of B, S2 becomes the purely spatial unit vector §2 which points along the 
spin direction. In the present calculation the results will be organized to explicitly exhibit the 
angular correlation §2 ■ 1. 

Energy and momentum conservation determines the allowed kinematical region in the variables 
E and Ei for process ©• This region, which is referred to as the Dalitz plot and is represented by 
the shadowed area depicted in Fig. [1] and labeled as I, is bounded in Ei by 

Ef" <Ei< ^r", (6) 

where 

pmax,min ^ {M2 + E T if + Mf 

^ 2{M2 + EtI) ' ^^ 

while the charged lepton energy falls within the interval 

m<E<Em, (8) 

where 




FIG. 1 : Kinematical region as a function of E and Ei for baryon semileptonic decays. The areas I and I + 1 1 
correspond to the Dalitz plots of the processes A^ B + I + Vi and A^B + l + Vf + '^, respectively. 



Similarly, area II in Fig.[T]is bounded by 



m<E<E^ 



where 



Er 



(Ml - M2f + m2 



(10) 



(11) 



2 (Ml - M2) 

The distinction between these two areas has important physical implications that should be 
clarified. Finding an event with energies E and Ei in area 1 1 demands the existence of a fourth 
particle which in our case will be a photon and will carry away finite energy and momentum. In 
contrast, in area I this photon may or may not do so. In consequence, area 1 1 is exclusively a four- 
body region whereas area I is both a three- and a four-body region. We will refer loosely to areas I 
and 1 1 as the three- and four-body regions (TBR and FBR) of the Dalitz plot, respectively. 



III. VIRTUAL RADIATIVE CORRECTIONS 



The method to calculate the virtual RC to the Dalitz plot of unpolarized and polarized decaying 
baryons has been discussed in detail in Refs. pi] and [141]. It can be readily adapted to our case here 
of nonzero §2, so only a few salient facts will be repeated now. The virtual RC can be separated 
into a model-independent part M„ which is finite and calculable and into a model-dependent one 
which contains the effects of the strong interactions and the intermediate vector boson. To order 
(a/vr) (g/Mi)°, the latter amounts to two constants {a/Tx)c and (a/7r)(i which can be absorbed into 
/i and Qi of Mq, respectively, through the definition of effective form factors, hereafter referred to 



as f[ and g[. Thus, the decay amplitude My with virtual RC is given by 



My = Mo + M„, 



(12) 



where 



and 



M, 



a 
Yti L 



Mo(/)+ M 



pi" 



M 



pi 



\mMiJ y/2 



C 

— = [ubWxUa] [uefiiOxVu] ■ 



(13) 



(14) 



The prime on Mq in Eq. (fT2l) will be used as a reminder that the effective form factors appear 
explicitly in this amplitude. Also, to order (g/Mi)° the amplitudes Mq and Mp^ in Eq. (fT3l) [but 
not in Eq. (fT2l) 1 are limited to contain only the leading form factors /i and gi. The calculation of 
the model-independent functions (p^E) and 0'(-E) shows that they formally retain the same form 
given in previous work ^. The hats over them denote they are now given in the center-of-mass 
frame of the emitted baryon B. These functions read 

'1 



0(E) 



P 



tanh"^ (3-1 



In 



m 



+ < 



— + - In ^ 

/3 2 m 



NDB) 



2/? 
1 + /?, 



+ -tanh-/3-- 



(15) 



3, Ml 

- In 

y 2 m 



(CDB) 



and 



^'(E) 



"-i 



tanh /5, 



(16) 



where [3 = l/E, L is the Spence function, A is the infrared-divergent cutoff and with CDB and 
NDB we distinguish the results for the charged and neutral decaying baryon cases. The divergent 
term in Eq. (fT5l) will be canceled by its counterpart in the bremsstrahlung contribution. 

At this point we can construct the Dalitz plot with virtual RC by leaving the energies E and Ei 
as the relevant variables in the differential decay rate for process ©. After making the replacement 
(0]) in (fT2l) . squaring it, averaging over initial spins, summing over final spin states, and rearranging 
terms we can express the differential decay rate as ^ 



dr 



V 



a 



dQ\A'Q + -{A[<p + A'l<l)')-S2-\ 



A" 



a 



TX 



(i'a^ + i'aV 



(17) 



' We shall also use hats over other expressions to emphasize that the center-of-mass frame of B is being used. 



where 



dn = - 

2 



2M2. (18) 



M2Y Gl dEdEidVtd^i 
Ml 



2 (27r)5 

Let us notice that this expression of dVt has some differences with respect to the one of previous 
work [|2|]. These are the factor 1/2, which results from averaging over the spin of the initial baryon, 
and the factor {M2/M1Y, which arises out of the Lorentz transformation to the new reference 
frame. To recover the unpolarized decay rate one makes the factor 1/2 disappear by inserting in 
Eq. (I2I) the operator S(— S2) = (1 + 75 /2)/2 instead of © and adding the result to (flTI) . 
The functions Aq and Aq, which emerge in the uncorrected amplitude Mq, read 

A'^ = EElQ^ - Ep^{pr - lyo)Q2 - lil - vm)Qz - Elvxh^Q^ + vlh^^Vi - lyo)Q^, (19) 

and 

i[; = ^pii/oQe + ElQ-r, (20) 

where y^ is defined as the scalar product pi ■ 1 and can be expressed as 

Pl + l'- Ef .... 

y^ = — 2^1 — ' ^^'^ 

and also, by energy conservation, the neutrino energy E^ is given by 

E^^ = E^-E- M2. (22) 

The Qi are new functions of the form factors and are listed in Appendix El The hat is used to 
avoid confusing them with the ones of Ref. M\. The functions A[, A'(, A'2, and Ag that emerge in 
these virtual RC read 

A[ = D^EEl - D2l{l - piz/o), (23) 

A'l = ^EElD,, (24) 

i'2 = E{1 - piyo)D, - lElD,, (25) 

and 

EMo 
A'^ = ^{l-Piyo)Ds, (26) 

where the coefficients Di are quadratic functions of the effective form factors. Explicitly, they are 

Di = fC + H\ D2 = fC - g[\ (27a) 

D, = 2{-g[' + f[g[), D, = 2{g[' + f[g[). (27b) 



Here we use also the effective form factors, so that our result is uniformly expressed. This is a 
rearrangement of second order ma/ix, which we are free to make within our approximations. 

To stress the parallelism with our previous work we have used the same notation, but there 
should arise no confusion. The expressions given here apply to the present case only. 

IV. BREMSSTRAHLUNG RADIATIVE CORRECTIONS 



To obtain the bremsstrahlung RC we have to consider the four-body decay 



A- ^ 5° + r + z7^ + 7, 



(28) 



where 7 represents a massive photon with four- momentum k = (u;, k) and fc^ = A^. To 
obtain the bremsstrahlung RC in a model-independent way we shall use the Low theorem ISi la], 
which asserts that the radiative amplitudes of order l/k and (k)'^ can be determined in terms 
of the nonradiative amplitude without further structure dependence. Then, we can express the 
bremsstrahlung amplitude M^ as 



Mb = Mbi + Mb2 + Mb3, 



(29) 



with 



and 



M 



Bl 



eMo 



21 -e 



+ 



2pi-e 



2l-k + X^ A2-2pi-fc 



M 



B2 






(30) 



(31) 



V2 "^^ " "2l-k + X^ 

Mbi contains terms of order l/k and Mb2 contains terms of order (/c)". Although M^s contains 
also some terms of order (k)^, it can be ignored because its contribution to the decay rate is of 
order q/Mi ||2|] so we do not need its explicit form. The infrared-divergent terms are all contained 
in Mbi- 

Next, we have to replace Eq. & in Eq. (l29l) . square the resulting M^, average over the initial 
spins and sum over the final spins and over the photon polarization. To perform the latter sum, 
we proceed in two ways. First, we can use the rule of Coester ||8|] to account for the longitudinal 
degree of polarization of the photon, namely. 



E/ N/ ,N u (a-k)(b-k) 
(e.a)(e-6) =a-b — 



(32) 



where cu'^ = k'^ + X^ (here k is the magnitude of k) and a = (ao, a) and b = {bo, b) are 
arbitrary 4-vectors. Second, the infrared-convergent contributions can be calculated using the 
usual summation over the photon polarization, namely, J2e(^ ■ a)(e ■ b) = —a ■ b, with to = k. 

After a standard calculation, the bremsstrahlung differential decay rate dTs corresponding to 
the TBR becomes 



dTn = dT'u 



rfrg), 



(33) 



(s) 



where dV'^ denotes half of the unpolarized decay rate whereas dT^ contains the spin of the 
emitted baryon. Explicitly, using the Di of Eqs. (|27al) and (|27bl) ^, they are. 



dT' 



B 



1 a G\ 2M2 

2^ 2 (27r)6 



M2 
Ml 



d^'Px dH S'k d^p, 



u ^4 



X 



+ 



El E uj El, 



S {pi-P2-l-Pu-k) 



{uo — pkx) 
1 



E(uj - (3kx) 



DiE„ {uj + 2E 



rn?uj 



Eiuj - I3kx) 



E(u - (3kx) 



u 



+D2 



k ■ p J 1 - 



m 



E{uj — (3kx) UJ 






(34) 



and 



dV 



B 



1 a Gl 2M2 
'2n 2 (27r)6 



M2 
Ml 



d^Pi d^l d^kd^p 

El E UJ El 



'S\pi 



X 



+ 



(3^1 - k^x^/uj^) 



{uj — pkxy 
1 



{D3EI ■ p, + D^EJ) 



E{uj - (3kx) 



D3l-Pu[uj + 2E 



m?uj 



P2-I -Pu- k) (S2 ■ 1) 



E{uj - (3kx) 



E{uj - (3kx) 



UJ 



+D4E, 



k-1 1- 



<s) 



m 



E{uj - (3kx] 



+ 



E 



UJ 



+ 1 



(35) 



The above expression of dT^ is specialized to the angular correlation §2 ■ 1. To achieve this we 
have used the replacement §2 ■ p — > (§2 • 1)(1 ■ p), with p = pi, k, pj,, which is valid over the 
Dalitz plot after all other variables are integrated. Both expressions (|34|) and (l35l) contain the 
infrared divergence in their first summand within the curly brackets, which we analyze in detail 
in Appendix IHl The A^ that appears in Eq. (|30l ) contributes in Eqs. (|34l ) and (|35l ) with linear and 
higher powers, which will become zero in the X^ ^ limit. 



^ For the sake of the uniformity of our results, the effective form factors /{ and g[ may be used here, too. Again, this 
amounts to a rearrangement of order (o/tt)^, valid within our approximations. 



The integration over the neutrino variables with (5^(pi — 1 — p,^ — k) in Eqs. (l34l) and (1351) is 
trivial, but leaves a nontrivial argument inside the last 5{Ei — M2 — E — E^ — cu), which allows 
one to integrate over the photon momentum. Without further ado, the resulting expressions can be 
cast into 



dv' = -dn 

^ TX 



A'Jo + (pi + p[)D, + (p2 + P2)D2 



and 



dT^^> = -dn S2 ■ 1 



A'Jo + (P3 + P3)^3 + (P4 + p1) A 



(36) 



(37) 



where dd, A[, Ag, Di, D2, -D3, and D4 were already defined. The integrals over the angular 
variables of the photon, x = k ■ I and (^k, and over y = pi ■ I are left to be performed numerically, 
except in Iq which contains the infrared divergence and is obtained analytically. Its result is found 
in Appendix El The functions pi, p[, p2, P2, ps, P3, P4, and P4 thus read 



(3WE n\ [\ r^'difk i-x^ 



p'l 
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P'2 
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dy dx dip 
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dy dx , 
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l-(3x 
2^ dpk 1 - x^ 



l-^x ^ E 
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dy i dx 1 dpk 
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where 



with 



D l-px 

piKv - yo) 



P + x{ 1 



UJ 



1-px E 



UJ 



D 



D = El- (pi - 1) • k. 



(38) 

(39) 

(40) 

(41) 

(42) 
(43) 

(44) 

(45) 
(46) 
(47) 

(48) 
(49) 
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We have obtained the bremsstrahlung RC to the differential decay rate to order («/7r)(g/Mi)°. 
In the next section we present the total differential decay rate by gathering both virtual and 
bremsstrahlung contributions together. 

V. FINAL RESULTS AND NUMERICAL FORM OF THE RADIATIVE CORRECTIONS 

We have reached our first goal: To obtain the complete RC to the Dalitz plot with the §2 ■ 1 
correlation to order (a;/7r)(g/Mi)° restricted to the TBR. The final result is obtained by summing 
the virtual RC, Eq. (fTT] ). and the bremsstrahlung RC, Eq. (|33l) . The latter is obtained when 
Eqs. (|36l) and (|37l) are put together. Thus 

dT = dTv + dTs. (50) 

We can rearrange this expression into a simple form as 

dV = dn \A'o + -Qj - §2 ■ i U'o' + -Qn] I , (51) 

In L vr J J 



with 



and 



0, = i;(0 + Jo) + i;'0' + (pi + p[)D^ + (p2 + p'2)D2, (52) 



e,, = i;,(0 + Jo) + i'2'0' + (P3 + P3)Ds + (p4 + p',)D,, (53) 

where all the ingredients are given in the previous sections. 

Notice that the only difference between the NDB and the CDB cases is found in the function 
(j){E), Eq. (fT5l) . The bremsstrahlung RC does not make any difference between both cases because 
of the order of approximation in this work. 

We now come to our second goal in this paper. This Eq. (|5T1) has triple integrals over some 
angular variables ready to be performed numerically. It requires that the numerical integrals in 
the RC be calculated within a Monte Carlo simulation every time E, Ei, and the form factors are 
varied, a task that represents a non-negligible computer effort. We shall now discuss a second 
form of the RC that should be more practical to use. 

For fixed values of E and Ei, Eqs. (|52|) and (l53l) take the form 

©r = a-rfl + brfl9l + Cr-gf, (54) 
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because they are quadratic in the form factors. The subindex r takes the values r = I, II. The 
second form of RC we propose consists of calculating arrays of the a^, br, and Cr coefficients 
determined at fixed values of {E, Ei) and that these pairs of {E, Ei) cover a lattice of points on 
the Dalitz plot. 

To calculate the coefficients a^, h,. and c^ it is not necessary to rearrange our final results to 
take the form (|54l) . One can calculate them following a systematic procedure. One chooses fixed 
{E, El) points. Then one fixes /i = 1 and gi = 0, and obtains a^, one repeats this calculation for 
gi = 1, /i = 0, to obtain c^. Next, one repeats the calculation with /i = 1, gi = 1, and from 
these results one subtracts a^ and c^, this way one obtains the coefficient bj.. The arrays of these 
coefficients should be fed into the Monte Carlo simulation. Within this simulation the repetitive 
triple integrations are reduced into a form of matrix multiplication. 

We may close this section by stressing that none of the forms of our RC results is compromised 
to fixing from the outset values for the form factors when such RC are applied in a Monte Carlo 
simulation. 

VI. DISCUSSIONS 

Our final result for the RC to the Dalitz plot of BSD with the angular correlation §2 • 1 between 
the polarization of the emitted baryon and the direction of the charged lepton is given in Eq. (|5TI) . 
It is valid to order (a/7r)(q/Mi)^ and it covers the TBR of this plot. It meets the requirements 
discussed in the introductory section, namely, it is model independent, it can be used in all charged 
assignments in different BSD, the charged lepton may be e^, /i^, or r^, and it is not compromised 
to fixed values of the form factors of the uncorrected decay amplitude. The finite terms that 
accompany the infrared divergence in the bremsstrahlung RC are given in analytical form. The 
other terms in this correction are presented with triple integrations ready to be performed. 

This result may be used in a Monte Carlo simulation of an experimental analysis. However, 
performing the triple integrations every time E, Ei, and the form factors are varied may represent 
a very heavy computer effort. A more practical use of our result is through Eq. (|54|) . Numerical 
arrays of the RC to the coefficients of the quadratic products of form factors may be first obtained 
at a lattice of points {E, Ei) covering the Dalitz plot and afterwards be fed in the Monte Carlo 
simulation. The computer effort within it would then be reduced to a sort of matrix multiplication. 

The procedure of this paper may be followed in the future to extend the calculation of RC to 
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BSD with the observation of polarization of the emitted baryon to cover the angular correlation 
§2 ■ pi and the FBR. The precision of RC may be improved, while still preserving their model 
independence, by including terms of order {a/n)(q/Mi)"' with n = 1. It may be the case that 
experimental analyses should be limited to the center-of-mass frame of the decaying baryon A. 
Our results should be then adapted to this frame. Each one of these possibilities requires further 
serious efforts. They should be attempted as the need for them arises. 
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APPENDIX A: THE Qi COEFFICIENTS 



The Qi factors contained in Eqs. (fT9l) and (|20l) are quadratic functions of the form factors. For 
the spin-independent contribution they read 



M2 
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El + Ml ^2 , ^1 - Miq2 



Mo 



-F^ + 



Mo 






Qi = vr(^? + G? - 2fiGi + FiF, - G1G2) + B1Q5, 



M^Qs 



2Mi 

mT 



(F1F2 + G1G2) + 



E, + Ml 2 ^ ^i-Mi 2 



Mo 



Mo 



2) 



(A3) 

(A4) 
(A5) 



13 



whereas for the spin-dependent contribution they read 



Qe 



E^-M^-(3piyo 



M2 

'El - PpiVo 

M2 

l + (l-/3') 



Fl + 



El + Ml - (3piyo 



2F1G1 + 



Mo 



l + {l-(3^ 



.2 (^PiVo , 



E 



Gi — — (-F1F2 + G1G2) 

iVl\ 
, E\ M2 + 2E- Ppivo 



E 



M, 



Ml 
E \ , M2 + 2E-ppiyo 



Ml 



F1G2 



Ml 



+ 



Ml 



F2G1 



+ {1-P' 



E 

'Mo 



l-^--]FiG, 
Ml Ml' ' ' 



M2 E , ^ 



+ 



^ + (^-^^)m2 + 



E_ El - /Spivo f _ 2{EiE - pilyo) + 2M2E1 
M2 \ Ml 



FoG 



T-J2 






E 
Ml 
E^ 
Ml 



Ml El E El- PpiVo 
M2 Ml Ml M2 
'El-^{piy,-l)- 



{F2G3 + F3G2) 



M, 



F3G3, 



(A6) 



and 



Qi 



1 + 



Ml 



El - Ml 

E 



Fl + 



1- 



Mi 
M> 



Ei + Mi 
E 



pI 



Gi + ^iFiF2 + GiG2) 



+ 



' ^ 'M2 E 



2FiGi H — F1G2 H — F2G1 



Ml 



Ml 



+ '^-^')i| 



^i-^ipr^ , Ei + Mi ^^- 

■riLrs H — r^Lri 



Ml Ml 

In these expressions we have used the form factors Fj and Gi, which read 



(A7) 



Mo 






F2 — —2/2, F3 — f2 + fs, 

G2 = -2g2, Gs = g2 + g3- 



APPENDIX B: EXTRACTION OF THE INFRARED DIVERGENCE 



Here we discuss the procedure we followed to identify and isolate the infrared divergence 
contained in the Iq function introduced in Eqs. (l36l) and (l37l) . We only show how the infrared- 
divergent term is calculated in the spin-independent contribution, because in the spin-dependent 
one the procedure is analogous. The term where the infrared divergence is contained is 



E |Mml 



e'G'y 



E 



spins, € 



21 -e 



2pi-e 



2l-k + \^ \^-2pi-k 



2M2 



Mimrrii, 



{N + N'), 



(Bl) 
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where the factors A^ and A^' are 

N = DiEEl - D2l{l - Piy), N' = -uo{EDi + IXD2), (B2) 

and 1/ = pi ■ 1 and x = k ■ 1. A^' is proportional to uj so that it is infrared-convergent, and it is 
absorbed into Eqs. (|38l) and (|40l ). Therefore, we here only consider the contribution of A^. With all 
generality, we can orient the coordinate axes so that the momentum of the charged lepton is along 
the 2;-axis and so that A is in the first or fourth quadrant of the plane (x, z). Thus, performing the 
sum over the photon polarization in the Coester representation and rearranging terms yields 

dn = -dn^ f dx r d^, /''" dk^-g{9,)N ^'^}-^lf'f\ (B3) 

TT 27r J_i Jo Jo ^ (^ - Pkxy 

where co"^ = k"^ + A^, km is the maximum value of the photon momentum. 



k ^°' - ^P^ - '^ (B4) 

and g{Oi) emerges from the last 5 function. It is given by 

,(90 = ^ ;'■"'' „ . (B5) 

a sm &i — cos Oi 

61 is the polar angle of the decaying baryon, a = 2pi(l + k cos 9k), b = 2pik sin 9k cos (fk, 9k and 
ifk are the polar and the azimuthal angles of the photon, respectively. 

We find it convenient to consider the partition (0, A A;), (A A;, k^) of the integration interval 
(0, A;^) for dk, with A A; arbitrary. Thus 



dV% = -dQ^-^ / dx / du^k / dk—gi9i)N- 

TT 27r y_i Jn Jn UJ 



2 



+ -dQ—^-^ / dx / d^k / dk—g{9i)N 
vr 27r J_i Jq J^k ^ 



(to - Pkxf 



= dT%^'^ + dT%^'\ (B6) 

We must stress that our result does not depend on AA;. The divergence is now in the first integral 
and because AA; is arbitrary, we can make it slightly larger than A, i.e., AA; > A. Then we can 
approximate g{9i) in the first integral, by allowing A; ~ 0. Thus, 

9(0,, ^ ^. (B7) 

Also, in A^ we expand y in powers of k up to first order, 

yc^y, + f'k + 0{e), (B8) 
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where 



/' = — [El + (/ - p,y,)x - pi(l - ylf/\l - x^fl^ cosO^] , (B9) 

so that N takes on the form 

N ^A'^ + D2Pilfk. (BIO) 

The first integral in Eq. (|B6I) becomes 

dTf'> = -dn^ dx d^k dk—--{A[ + D2Pilf'k)- -^, (BU) 

and then the infrared divergence is finally contained in the first summand of the above equation. 
The second summand picks up a factor of k so it is infrared-convergent and we can use tu = A; in it. 
It will cancel away with one term of dF'^ [see comment after Eq. (|B21)) ] . A further simplification 
is obtained by identifying in Eq. (|B1H) the integral of Kinoshita and Sirlin ['9'] 



/ dx 

27r 2 



difk I dk- 
io 



UJ {uj — pkxy 



2 In 






/tanh"^/3 



C, (B12) 



where 



C' = 21n(/i-2) + l + -/i 



2 + In 



l-P' 



4^ 



2/3 • 
1 + /?, 



2 ^ 2 ' 



(B13) 



and 



Ji = - tanh"^ /3. 

P 



(B14) 



Let us now analyze the second summand in Eq. (IB6L (iF^^ '. We shall change the integral over 
k into an integral over y. For this purpose, we rewrite the photon momentum as 



2D' 



(B15) 



and 



where 



N = A[ + -D2F, 



(B16) 



F = 2p,l{y-yo), 



(B17) 



and, rearranging Eq. (|49l ). 



D = ^° + /a;-pi(xy + a/(1 -y2)(l -x2)cosv2fc). 



(B18) 
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Then we can replace 



dk 



dy 



(B19) 



2Dgie,y 

The integration limits change to y{Ak) ~ yo + f'^k and y{km) = 1- Notice that the upper limit 
km of (|B6I) is replaced in this variable hy y = 1. 
Putting all these changes together yields 



dvT = -dCi^^ 

^ TT 27r 



1 p2tt p1 

dx I dipk / dy 



1 







yo+f'Ak 



A\ 1 
— + -D2 
F 2 



X 



{i-(3xy 



(B20) 



and the integration over y gives 



dr 



ir(2) 
B 



—dCl-—A[ / dx— — -^ 



27r 



civ^fcln 



Mi(l-yo) 
MJ'Ak 



IT 2n 



2lT 



dx / dipk-D2 
1 Jo ^ 



[l-y,-fAk)- 



l-x' 



(B21) 



(1 - Pxf 

We notice that the term proportional to /'A A; in the second summand of this latter equation cancels 
precisely the second summand proportional to D2k of Eq. (IBllI) once the integration over k is 
performed in it. 

Finally, the resulting expression for dV^^ becomes 



dT 



B 



dTf'> + dVf^ 



-dnA[ <! 2 

+ -dQA[ 



a 



IT 



dM[ 





tanh ^(3 ^ 


i 


2 

taa 


(3 
'tanh"^ p 


I 


[ P 



p 



In 
-1 
In 



+ C 
Mi{l-yo) 



'Ak 

X 

In 

[ Ak 

Mi(l-yo) ' 
A 



+ Ci 



DoQ 



+ C + Ci 



2<-^2 



+ -dQD2C2 



-dh 

TT 



A[lo + D2C2 . 

The function Iq is defined in Eq. (|B22I) as 



(B22) 



'tanh"^ p 



-1 



In 



" Mi(l-7/o) 
A 



C + Ci, 



(B23) 
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where 
Ci = 



2tt 



P' r- r 1 

■— / dipk I dx 



X 



-1 



{l-Pxf 



-- i In 



(a+- l)(l-a-) 



In 



4pi/Mi 
(a+ + l)(a- + l) 



4pi/Mi 



In [Ml/'] 
/3(l-/3xo) /3 /3"^ 1 + /3 

H 7^ ^ - In :; h (1 - Xq) 



/9 



^-^m^ 



/?Xo 



-{l + xo) 



81n2 + ^-^ln(l + /3)- 



L/?(l-/?Xo) /5 p-l 

ln(l - /3) + 2 + 



/?Xo 



-^ln(l-/3) + 2 + ^-^-^(l-xo) 



- Xo) ln(l - xo) 



21n(l -/5xo) + 2 /^(^ + ^") (i + a;„)in(l+xo) + ^ 
1-/5x0 /5 



L 



1-/3 
1 -/3xo 



-2L 



1-/3 
1 + /3 



1 + /5 / 1 + /3 1 + /5 2 \ 1+13 



where 



Xq 



I - PlVO 



and 



E^uTPi 



(B24) 



(B25) 



and C is given in Eq. (|B13I) . It possesses the right coefficient to exactly cancel the infrared- 
divergent term in its counterpart in the virtual RC, Eq. (fTSl) . 
On the other hand, C2 given by 



C2 =Pil—{l -yo) 



dx- 



1-x^ 

i-pxy 



(B26) 



is absorbed in p2 and pa. 
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